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O [ Abstract 

J>y We extend the results of Ref. [JJ on one-photon electric dipole tran- 

sition line shift and broadening to the case of two-photon transitions. 
As an example we consider the laser induced transition in antiprotonic 
helium produced in helium gas target. The transition is between an- 

^. ■ tiprotonic helium states (n, I) = (33,32) and (31,30). 
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1 Introduction 

In the present paper we propose an approach to the evaluation of the density 
shift and broadening of the two-photon transition line profile. The foundation 
of our considerations is the method developed in Ref. [1] for ET-transitions. 
We slightly modify the mathematical manipulations used therein and derive 
close expressions for the shift of the resonance frequency and line broadening 
in terms of the perturbing potential mean value. 

Our goal is to reduce the experimental uncertainty of recent high-precision 
spectroscopy measurements [2J of antiprotonic helium QoHe), related to the 
effects of the collisions of antiprotons with helium atoms. The example which 
we consider thoroughly is the two-photon transition from the initial pHe state 
(n, /) = (33,32) to the final state (31,30) [3]. The transition is induced 
by external monochromatic electromagnetic waves and is influenced by the 
helium gas target. Unfortunately, the available data about the pHe - He 
interaction potential [1J do not cover the whole range of interparticle distances 
which we are interested in. This forces us to use an extrapolation of the 
known potentials. Two possible extrapolations are considered and compared. 

2 Two-photon transition in low density gas 

The general footing of our consideration is as follows: A quantum system 
(emitter) is subject to a perturbation due to its randomly propagating neigh- 
bors (perturbers) and also interacts with an external electromagnetic field. 
The time dependent Hamiltonian of the entire system is 

H{t) = H + V(t) + W(t) (1) 

where H Q is the unperturbed Hamiltonian of the emitter, V is the emitter- 
perturbers potential, and W is the electromagnetic interaction. 

We suppose that the emitter possesses a full system of discrete states 
{\j >} with energies {Ej}, i.e. H \j >= Ej\j >. In what follows we shall 
distinguish three of these states, namely the successive states \i >, \m > and 
|/ > such that Ei > E m > Ef. The state \i > will be our initial state and 
|/ > will be the final one. 

In our consideration the electromagnetic field is a superposition of two 
plane waves with frequencies Uk, k — 1,2. The emitter - electromagnetic 



field interaction W is of the form 

W(t) =vf k cos(u k t). (2) 

The coefficients w k are emitter depending. For example, if the emitter has a 
dipole moment d then w fc = d.E k where E k are the electric field amplitudes. 
It is supposed that the frequencies u k are different but close to the resonance 
ones 

Ui PS U mi = (Ei - E m )/h 

U 2 PS Uf m = (E m - E f )/h. (3) 

We make four assumptions for the potential V pp. We suppose that, first, 
V it is too weak to cause quantum excitation in both emitter and perturbers. 
Second, the target density is low enough, so that the emitter interacts only 
with one perturber at a time via the pairwise emitter - perturber interaction 
V°. Third, V° depends only on the distance R between the emitter and per- 
turber. Forth, we adopt the Anderson approach jl] and treat the perturbers 
classically. Moreover, following [lj, we assume that the perturber's trajectory 
R(t) is entirely determined by the potential V° (plus initial conditions). 

There are two important consequences of the above assumptions. First, 
the perturbing potential V is a sum of pairwise interactions 

V(t) = J2V°(R n (t)) (4) 

n 

where n runs over all perturbers. Second, both V° and V have the set {\j >} 
as an eigenvector system (with v°(i?) and Vj(t) = J2n V j(-Rn(^)) as eigenval- 
ues). Therefore, these operators commute with the Hamiltonian H and the 
evolution operator of the initial system plus perturbation is diagonal with 
respect to the system of states {\j >}. Its matrix elements are determined 
by the action over the corresponding classical solution: 

< k\u Ho+v (t f Mj >= V"* i * (t '"* ) ~ M *' A) - (5) 

Here 

Vj(tf,t i )) = l;£ f dTV j (t) (6) 

is the action of the perturbing potential when the emitter is in state \j >. 



The matrix element of the entire system evolution operator U(tf,ti) be- 
tween states \i > and < /| up to the second order is 

</|*7(t/,ti)|i> = e^ {t f- u)E f /h+iri f {t ^ u) [5 fi 

h Jti 

__ e iu fiU f tf dTl ^ e -^n-%i(nA) f/ .( ri ) 
f 1 dr 2 e-^^-^^'^Wjifo)]. (7) 
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Here tjm = f]i — f]k and Wki is the matrix element of the electromagnetic 
interaction operator W between states < k\ and \l >. We can simplify 
eq.([7|). We can neglect the common phase factor. When i ^ f we can omit 
the diagonal term as well. Keeping on only slowly oscillating terms we can 
reduce the sum over the full system of intermediate states to one term only, 
namely \m >< m\. Finally, we can neglect by the same reason the first order 
term. As a result the evolution operator matrix element between different 
states < f\ and \i > takes the form: 

< f\U{t f ,ti)\i > = W/m 7 mt f tf dn e >/»-*)«-^(iA) 

x / dri e~ i ( UJmi ~ u!l * >T2 ~ irimi ( T2 ' ti * > (8) 

Jti 

It is clear that the line shape, shift and broadening are determined by 
the double integral in eq.flH]) and do not depend on the factor Wj m w^/^. 2 . 
In what follows we shall neglect this factor and concentrate our attention on 
the integral which in the limit U — > — oo tf — > oo we denote by U 



U = dt e-^fm-^t-ivfUt)}^ ( 9 ) 

I(t) = f dre^ Wmi " wl)T "* w W. (10) 

A few remarks. What we really know is v®(R) -- the potential between 
the perturber and emitter in state |/ >. With proper boundary conditions 
(e.g., impact parameteiLj r and velocity v) we can, using v°(i?), calculate the 



Hhe distance between the emitter and the straight line which coincides with the per- 
turber's trajectory at infinity 



perturber's trajectory R(t). Knowing the trajectory we can find v®(t;r,v) 
but this is not enough. What we actually need is v/(t) which is a result 
of successive random pairwise interactions. To find it we have to choose 
a sequence of interaction moments {£«} so that V(t) = J2iV°(t — U;ri,Vi). 
Using V(t) we have to calculate the probability of the process of interest 
and then to average over all possible time sequences, impact parameters and 
initial velocities. The task seems hopelessly complicated but before to make 
some simplifications let us see what we can get. We expect, on the base of our 
assumptions about V(t), that both r)f m (t) and r] mi (t) are stair-like functions 
of time, i.e. (suppressing for a while the state indexes of rj) 

r](t)=r] + c t + fi(t). (11) 

The important term here is cot while fj(t) gives a small irregular variation 
around it. (The constant i]q is irrelevant. It gives an overall phase in the 
transition amplitude.) Let us introduce the following simplifying notations: 

U = Uf m - w 2 + (c ) m - (c )/ (12) 

u = u mi -uj 1 + (co)i - (co) m (13) 

V = Vm-Vf (14) 

V = Vi-Vm (15) 

where (cq)j is the cq coefficients oirjj. We will need also the Fourier transform 
of fj and fj 

#) = |^e^e(o, m = Jdt^m)- (is) 

After regrouping the leading linear in t terms in I(t), the integral takes the 
form 

/(£)= (' dre- iWT - if,{T) . (17) 

Now, because V is small we suppose that fj(t) and fj are small as well (see 
the comment below about this point), so we can make series expansion of 
e~ tf> and e~ 1 ^. Therefore, I(t) and U can be rewritten as follows 

I(t) = f dr e- WT (1 - ifj(r)) 

-iut ( 1 , f A c e (0 Jit 



* \z, + }«^r) (18) 



U = Jdt e -^+-)* (l - ifjit)) f 1 + Jdi ^zje^ 
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= Jdt e-^* (L + / d£ (2& + Ml) e*) 

\ w a; a; / 

The line shift A of the two-photon transition can be read immediately from 
eq. (fT9l) . Recalling the definitions of u and u> we see that u + u = (Ei — 
Ef)/% — co 1 — co 2 + (c )i — (c )/ and therefore 

A = (co)i - (co)/. (20) 

Note that 

(co)i= fe4// rv ' (r) ' (2i) 

i.e., Co is the mean value of the perturbing potential. 

It is easy to find the line shift, but the situation with the line broadening is 
more complicated. At the present moment we can not say anything about it, 
because we do not know the explicit form of 0(0 and 0(£). We shall devote 
the rest of the paper to show that for small but non-zero £ the functional 
form of (and G) is entirely determined by the coefficient Co and it is 

6(0 oc |. (22) 

We want to step back a little and comment the expansion of e~^^ used 
in eqs.( H8fE?|) . The assumption that rj is small for any t is correct, if we 
have a gap (or cut off) near the zero in the spectrum of fj. The easiest may 
to ensure the existence of a gap is to suppose that V(t) is periodic. The 
idea is that the mean value of the perturbing potential which determines 
both the shift and broadening of the line is independent, according to the 
Central Limit Theorem, of the assumption for periodicity, but can be easily 
calculated using it. 

Supposing that the perturbation is periodic its period to is the time be- 
tween two emitter - perturber impacts. This time can be determined by the 
density N and temperature T of the perturbers 

*o = AT 1 - ( 23 ) 



where v = J8kT/irm is the mean speed of the perturbers (m is the reduced 
mass of the system) and r max is the distance above which we can neglect the 
perturbation. Obviously, r max depends on V° but is also in our hands. 

The assumption of periodicity leads to a dramatic simplification in the 
calculations. In the periodic picture 

V{t) = y £V°{t + nt ) (24) 

n 

where Vo(t) is the mean emitter - perturber potential 
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v{t) = nz dr2r tAw 



x dve- mv2/2kT v 2 V°(t;r,v). (25) 



Certainly, Vo(t) commutes with the Hamiltonian. We denote its eigenvalues 
with v°. The time parameter in eq. (|25j) is chosen so that t = corresponds 
to the apex of the perturber's trajectory. Therefore, vf(t) and v/(i) are 
symmetric functions of t. For Vi(t) we can write down a Fourier series 

vi(t)/n = (c ) i + 2(c fc )iCos(A;u;ot). (26) 

1 rto/2 
^)i = 7T / dt cos(nu t)^(t). (27) 

toll J-to/2 



Note that, because of the existence of r max , what we really calculate is 

rfi/2 

i/2 



^)i = 7t/ 1 dt cos(kuj t)v°(t) (28) 

t n J-ti/2 



where t\ is the time for which the perturber propagates through the area of 
nonzero potential. According to one of our assumptions listed above to >> t\. 
Eq. (128]) is an origin of a very useful symmetry 

to - y o&q 
V°(t) -)• aV°(t) (29) 

provided ato > t\. Proof: Let us change to so that tQ ew = t /k. Therefore, 
to™™ = kuj and 



— - / cos(nkuot)v° (t) .dt = kctn (30) 

ton J-ti/2 
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Coefficients c™ 6 "' and c^ n correspond to one and the same frequency nuQ SW = 
kncjQ which means that c new (u>) = kc(u>). If we rescale both i and V° as it 
is prescribed by eqs. (l29j) then we will get that c new (u) = c(u). Now look 
at eqs.([23~[ [251) . We see that different choice of r max leads exactly to the 
transformation ([29]) . Increasing r max we get weaker mean perturbation but 
it happens more often with the same gross effect. When 1/a is not integer 
there will be a shift in the overtone positions, but still the Fourier coefficients 
will lay on one and the same curve. This curve - the envelope of the Fourier 
coefficients is the important one for us and it determines the line shape. Note 
that symmetry ([29]) holds both for a < 1 and a > 1. In the latter case the 
only limit on a is determine by the condition V° /h << 1 no matter to 
what r max it corresponds. We shall use such transformation to probe the line 
shape at small frequencies. 

Using eq. ( [26]) we get the following expressions for the actions r\ and r\ (see 
also eq.flnj): 

rj{t) = cot+Y.P^M^ot)) (31) 

V(t) = c t+J2l^Mkoo t)). (32) 

In the right hand side of the above equations all therms are smalu but the 
first ones (because t is arbitrary). Therefore 

e -iv(t) = g-fcot^ _ t pL S i n (kuj t)) (33) 

kOOq 

and 

-' ^ (J-V.) 2 ) sin(fa ° f> ) (34) 

Now we shall need some well known formulas 

"° dt e~ iuJt cos(wt) = Tr(5(u + cu) + 5(uj-co)) (35) 

oo 
oo 

dt e~ lult sin(wt) = in {5{u + u) - 5{u - u)) (36) 



2 this gives another way to determine the maximal a (or the minimal u>o) we can use 



Using these equations we obtain that U (up to a phase) is 

U = 2ir(-6(u + u) + 

\uj 

+ ~r^\ — v — r^ — v )6(u + u- ku ) + 

+ f /" - ( ~\ s ) <% + ^ + hjo)) (37) 

Hereafter we suppose that one of the frequencies u k , say a; 1 , is fixed. Then the 
probability interpretation of eq. ( 1371) is exactly the same as for the simple one- 
photon transition amplitude between unperturbed states. The only difference 
is that instead of one line now we have a bunch of closely separated lines 
with intensities proportional to the square of the coefficients in front the 
delta functions. In practice, where the real potential is not periodic, we see 
the envelope of these lines. This is in agreement with eqs. (128I29D according 
to which at the limit to — > °° the inter line distance Uq is so small that the 
Fourier coefficients c^ form a line. The important moment is that the Fourier 
coefficients calculated for any t lay o n this line. Therefore, we can fix to, find 
the Fourier coefficients {ck} and {c k } and interpolate them with functions 
c(f) and c(£). 

The expression for U given in eq.( l37j) is very close to that in eq.(IT9~l). It 
is clear that the line shift again is determined by eq.( l20p but now we can 
say something more about line shape. The envelope functions c(£) and c(£), 
according to eq.()28l). behave like constants for £—*•(). So, we can smoothly 
continue them for negative £ making them symmetric with respect to £ = 0. 
As a result U for sufficiently small but non zero £ is (see also eq.( 



Uoc^—f^-^-) (38) 

(u+u) \U CO J 

Eq.( l38"]) allows us to estimate for given u 1 the line broadening a (calculated 
as the difference between frequencies for which the probability is half of its 
maximum) 



a = \ [u 



±V2c ) 2 ±AV2c u. (39) 

The choice of the sign in eq. (139]) depends on which of discriminants is pos- 
itive. If both discriminants are positive then we have double line. (The 



doublet can be seen easily as one much broader line in experiments with non 
monochromatic light sources.) 

The dependence of the Fourier coefficients c on density N can be read 
from eqs. fl23ll28|) 

<*W = (^) »■ («) 

As a consequence of eq. (1401) the line shift and broadening also depend on N 

A(A0 = (5ffl + flffl)jv (41) 



N, 



o 



' V2\co{N ) 
a(N) \u>\ if |c |, |c | << \u\. (43) 



cr(N) = ^ U; ' ]N if Ifiol, |c | » |w| (42) 



3 Numerical results 

The system we consider [5] consists of antiprotonic helium as emitter and the 
helium atoms in a gas target as perturbers. The initial, final and intermediate 
states are \i >= (33,32), \f >= (31,30) and \m >= (32,31) respectively. 
The target is at p — 1 mbar and T = 6°K. 

We use two sets of perturbing potentials {v°, v^ and v° } corresponding to 
two different extrapolations of the data we have about the potential energy 
surface (PES) for the pHe - He interaction. The results obtained from the 
first set of potentials will be indicated by prime and those obtained from the 
second set by double prime. 

The double integral in Eq. ( |25|) was calculated as a left Riemann sum over 
a regular set of iV = 1000 points for the impact parameter r < r max = 25 
a.u. and using a Gauss-type quadrature formula with M = 6 points for the 
average over the Maxwell distribution for the velocity. 

We obtain the following values for the coefficients cq and c needed to 
estimate the line shift and broadening according to eqs. (j4"Tl - |4"3"|) 

c' (N )/N = 5.2 10- 13 [Hzcm 3 ] (44) 

do(N )/N = 4.8 10- 13 [Hzcm 3 ] (45) 

c'o(N )/N = 1.4 10" 12 [Hzcm 3 ] (46) 

c£(N )/N = 1.4 10- 12 [Hz cm 3 ]. (47) 



This gives that the line shift is 

A'(iV )/iVo = 1.0 l(T 12 [Hz cm 3 ] (48) 

A"(N Q )/N Q = 2.8 l(T 12 [Hz cm 3 ]. (49) 

We can use eq. (j4"2"|) to obtain the line broadening when |o;i| « |u; m j + Co|, i.e., 
when we have a fine tuning between the first laser frequency and the inter 
level distance between the initial and intermediate states. Then the result is 

a'(N )/N = 6.8 l(T 13 [Hzcm 3 ] (50) 

a"(N )/N = 2.0 10~ 12 [Hz cm 3 ]. (51) 

Eq. ( I43p describes the off-resonance situation. In this case the line broadening 
does not depend on the target density 

a(N) = \u mi -ui\ if \uj mi -wi| >> |c |, |c |. (52) 

In every other case the general formula (1391) for the line broadening has to 
be used. 

The discrepancy between results for the two approximating sets of poten- 
tials is about a factor of three. This indicates that the PES extrapolation we 
have used in the construction of the potentials is not reliable. The problem 
could be solved only by extending the PES to shorter distances between the 
antiproton and He + ion, corresponding to the average radius of the antipro- 
ton orbit in pHe in states with n ~ 30. 
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